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^ Abstract 

04 , We study the decorrelation in azimuthal angle of Mueller-Navelet jets 

'n— «' ' at hadron colliders within the BFKL formalism. We introduce NLO terms 

in the evolution kernel and present a coUinearly-improved version of it for 
all conformal spins. We show how this further resummation has good 
' convergence properties and is closer to the Tevatron data than a simple 

LO treatment. However, we are still far from a good fit. We offer estimates 
JL ■ of these decorrelations for larger rapidity differences which should favor 

JjJ the onset of BFKL effects and encourage experimental studies of this 
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observable at the LHC. 



1 Introduction 



X 

^ In this paper we continue the analytic study of the azimuthal decorrelations in 

Mueller-Navelet jets using the Balistky-Fadin-Kuraev-Lipatov (BFKL) equa- 
tion [I] beyond the leading order approximation initiated in Ref. |2j. We in- 
vestigate the inclusive hadroproduction of a pair of jets with large and similar 
transverse momenta produced at a large relative rapidity separation Y [3]. In 
principle, when this rapidity Y between the most forward and most backward 
jets is small there is not much phase space for the production of extra radiation 
and a fixed order perturbative calculation should be enough to describe the ob- 
servable. However, if Y is large enough to make the product OgY ~ 1 then a 
BFKL resummation of these terms to all orders should improve the accuracy of 
the calculation. 

In Section [2] we define our observables at partonic level and argue that the 
particular choice of rapidity variable removes the dependence on the parton 
distribution functions in normalized cross sections. The explicit expressions for 
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angular differential cross sections including the next-to-leading (NLO) BFKL 
kernel |4j with full angular dependence are also presented. In Section [3] we 
study in detail the structure of the scale invariant NLO BFKL kernel for dif- 
ferent conformal spins. We show how the convergence of this kernel in terms of 
asymptotic intercepts is poor for zero conformal spin while being much better 
for larger ones. After this, a prescription to improve the coUinear structure of 
the kernel using a shifted anomalous dimension for the full angular dependence 
is derived. The main consequence of this study is that even though coUinear 
poles are removed for all conformal spins only the asymptotic intercept corre- 
sponding to the angular averaged case is modified, while the intercepts of the 
angular dependent components are hardly affected by the resummation. In Sec- 
tion [4] we use these kernels to make predictions for hadron colliders. We further 
discuss the need of a coUinear resummation to generate stable results against a 
change of renormalization scheme and show that the resummed results provide 
a better description of the Tevatron data than those using LO BFKL. However, 
our results provide too much decorrelation when they are compared to the ex- 
perimental measurements. From the theoretical side this could be due to the 
approximations needed to obtain our analytic results, e.g., the choice of rapid- 
ity variable and use of leading order jet vertices. In principle, our predictions 
should be more reliable at larger rapidities where multiple emissions between 
the two tagged jets are favored. It would be very interesting to further study 
these Mueller-Navelet jets at the Large Hadron Collider (LHC) at CERN to 
gauge the importance of configurations in multi-Regge and quasi-multi-Regge 
kinematics and try to investigate other possible jet topologies dominated by 
them. This point is further discussed at the end of our work, in the Conclusions 
section. 

2 Cross sections 

We are interested in the study of normalized differential cross sections which 
turn out to be quite insensitive to the parton distribution functions. Therefore, 
to a good accuracy, the present analysis can be performed at partonic level 
and we focus on the parton -I- parton jet -I- jet -I- soft emission process. The 
differential cross section is then 



where as = asNc/ir is the strong coupling, gi.2 are the transverse momenta of 
the tagged jets, and Y their relative rapidity. The influence of the distribution 
functions is larger if the exact definition of the rapidity difference as in Fig. [1] 
is considered. For convenience we take Y as a fixed parameter since this allows 
us to perform all the necessary Mellin transforms and proceed with our analysis 
analytically. In more detail, the choice Y = hi{xiX2s/so) corresponds to a 
change of the energy scale sq — ^gfgf to a fixed value. This change can be 
also understood as a NLO contribution to the vertices coupling the tagged jets 
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to the external hadrons. The uncertainty associated to the choice of sq will be 
taken into account in the phenomenological discussion at the end of Section |4l 




Figure 1: Representation of Mueller-Navelet jets at a hadron collider. 
The gluon Green's function 

/ (91 , 92 , Y) = J ^e-Y/. (gi ,92), (2) 

carries the full dependence on Y and fulfills the NLO BFKL equation which, in 
the transverse momenta operator representation 

qWd = QtWi) , (3) 

with normalization 

(gilikl) = (9*1-92), (4) 

can be written as 

(uj - a,Ko - alki^ U = i- (5) 
This representation is useful when it acts on the basis 

(9>,n) = (6) 

which includes the dependence on the modulus and azimuthal angle of the differ- 
ent emissions. The Mellin-conjugate variable of the modulus is the real parame- 
ter ly, and the Fourier-conjugate parameter of the angle is the integer conformal 
spin n. 

As the rapidity difference increases the azimuthal angle dependence is mainly 
driven by the kernel. This is the reason why, in the present work, we make use 
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of the LO jet vertices which are much simpler than the NLO ones calculated in 
Ref. [5j . Following Ref . [2] we can then write the differential cross section in the 
azimuthal angle 4> = 9i — 62 — tt, where 0i are the angles corresponding to the 
two tagged jets, as 



da (as,Y,pf 2) 



J2 e™^C„(Y), 



n— — 00 



with 



and 



C«(Y) = — 
Zn 



dv 



Pl\ pX{\n\,l+iv,as{piP2))Y 



f^Q Xo 7) 
8N, 7(1-7) 



(7) 



(8) 



(9) 



Throughout this work, the coefficients Cn are not evaluated at the saddle point, 
but obtained by a numerical integration over the full range of v. In the above 
expression the LO kernel, Kq, has as eigenvalue the function 



Xo 



(n,7) = 2^(1)-V(7+|) -V^(l-7+^) , 



(10) 



with tp being the logarithmic derivative of the Euler gamma function. The last 
term in Eq. ([9]) stems from the scale dependent part of the NLO kernel, i.e. from 
the running of the coupling. Its explicit form, in our representation, depends on 
the impact factors and is discussed in more detail in Refs. [Sill]- The action of 
the scale invariant sector of the NLO correction, Ki , in the MS renormalization 
scheme, on the basis in Eq. © explicitly reads [8] 



Xi ("•,7) 
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4sin^ (717) (1-27) 
1 ^ 



iy\ 2 + 37(1-7) 
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Nf) 2 (3 -27) (1 + 27)^ 



(11) 



where we have used the notation S = (4 — tt^ + bjSo/Nc) /12, /3o = (H-ZVc — 
2ny)/3 and C('^) = X^p^iP"" the Riemann zeta function. The (j){n,^) func- 
tion is of the form 



<^(n,7) 



E 



A: + 7 + 



fe=0 ■ ' ' 2 
+ (-l)('=+i) {13' {k + n + 1) + 13' {k + 1)) + 



V''(fc + n+l)-V/(fc + l) 

i/'(fc + 1) - i/'(fc + n- 



A: + 7+^ 
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with 



These scale invariant eigenvalues have a very interesting structure which we 
will study in the next section. Before this we would like to indicate that the 
full cross section corresponds to the integration over the azimuthal angle of the 
differential expression in Eq. ([7]) . This implies that it only depends on the n = 
component: 

a{as,Y,pl,) = ^lM^Co(Y). (14) 

In this paper we are interested in those distributions which are sensitive to the 
higher conformal spins. In particular, the average of the cosine of the azimuthal 
angle times an integer projects out the contribution from each of these angular 
components. It can be obtained using the ratio 

(cos(m0)) = (15) 



The associated ratios 



(cos(to0)) _ Cm(Y) ^^g^ 



(cos (n0)) C„ (Y) 

are also of interest since they can be used to remove the uncertainty associated 
to the hard pomeron intercept, i.e. the n = component, which we will analyze 
below. To study the behavior of all the angular components together it is useful 
to use the normalized differential cross section on the azimuthal angle: 



= — 1 1 + 2 cos {ncj)) (cos (n0)) 



- Y.T.J-'W) = ^^1 + O.cos(n^)(cos(n0))^ (17) 



Before comparing the results stemming from these expressions with the exper- 
imental data it is important to first analyze the convergence of the kernel at 
NLO for the different conformal spins. We proceed with this study in the com- 
ing section. 



3 Collinear resummation 

It is well-known that the BFKL resummation presents an instability when the 
NLO corrections are taken into account, for details regarding this point, see, 
e.g., Refs. [6l[9l[10]. For the observables studied in this paper we have found 
that if we use the NLO BFKL kernel as it stands the cross sections are very 
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dependent on the renormalization scheme. In particular, the term proportional 
to xo in Xi can be removed by a shift of the Landau pole of the form 



Agb = A^eT^ 



(18) 



This defines the so-called gluon-bremsstrahlung (GB) scheme [11] [T^ which 
is commonly used when dealing with soft gluon resummations. In this new 
scheme we find that some of our distributions even change sign and become 
unphysical. This is a manifestation of the poor convergence of the series. A 
crucial ingredient to improve the convergence of the perturbative expansion is to 
demand compatibility of the BFKL kernel with renormalization group evolution 
to all orders in the limit of deep inelastic scattering. This can be achieved if 
corrections to all orders are introduced by means of a shift in the anomalous 
dimension [13] . So far these types of resummations have been performed for a 
BFKL kernel averaged over the azimuthal angle and therefore they only affect 
the zero conformal spin sector. For our purposes in this work we must study 
the convergence of the eigenvalues of the kernel for all angular components. 

The renormalization group improved kernels are based on a proper treatment 
of the coUinear region of emissions. This region is insensitive to the azimuthal 
angle so we should not find a big effect beyond the n — case when we re- 
sum. We will see that this is the case because the asymptotic intercepts for 
the different angular components with n > are very stable under radiative 
corrections. 

To start our investigation it is necessary to extract the pole structure of the 
different contributions to the kernel around the ^ = —^,1 + ^ points. These 
poles are as follows: 
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(21) 
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Here II„ stands for the harmonic number [n + 1) ~ t}} [1) . 

The term due to running coupling effects — , in Eq. ^ , introduces 

the following modification of the single and double NLO poles of the original 
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kernel: 



j3o_ f 1-^° 
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(24) 
(25) 



There is sonie freedom in the way the resummation of colhnear terms can 
be performed. We find that the most natural scheme is an extension of that 
discussed in Ref. [TUj , which was first proposed in [B] , to include the dependence 
on all conformal spins. For completeness we have checked that the numerical 
results we will present are very similar for different resummations [?]• In this 
way, to obtain a convergent series for all values of the conformal spins, we use 
the prescription 



uj = as(l+^„a,)<{2V(l)-V(7+^ + ^ +'8„a, 



(26) 
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where the An and Bn coefficients are related to those of the original NLO kernel 

by 



An = CLn Ip' {n -\- 1) , 



(27) 
(28) 



It is worth noting that the solution to this transcendental equation can be 
approximated to a very good accuracy by the expression 



asXoiHn) + 0!l{ xi(|n|,7) 



00 Xoin,!) 
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This is an extension to the present case of the "All-poles" approximation de- 
veloped in Ref. JLQJ. 

In the presentation of our resummed kernels and experimental observables 
we will be using the MS renormalization scheme. We have checked that these 
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results do not significantly change when the GB renormalization scheme is used 
instead. This gives us confidence on the stability of our calculations. In Fig.[2]we 
have plotted the eigenvalue of the scale invariant sector of the BFKL kernel for 
different conformal spins. The value we have chosen for the coupling constant 
is Ois = 0.15. For simplicity, in these plots we have not included the term 
related to the LO impact factors. As a general feature we see how for non-zero 
conformal spins the LO and NLO kernels contain poles which are moving away 
from the < 7 < 1 region. These singularities are removed when the coUinear 
resummation, which is indicated as "Shift" in the plots, is introduced. We also 
show how the approximation provided by Eq. (|29p . denoted as "All-poles", is a 
very good one. A remarkable feature of the NLO kernel takes place for conformal 
spin 2. Here the terms proportional to 5'^^ in Eq. pT|) generate poles at 7 = 0, 1 
which we choose not to resum away. The structure for the higher n's is the same 
as for n = 1,3 with the poles "traveling" towards the left and right directions 
in 7. In these plots we can already see how the region relevant for asymptotic 
intercepts, around 7 ~ 0.5, is only sensitive to the coUinear resummation for the 
angular averaged component, n = 0. For the other contributions the intercepts 
are practically invariant under the introduction of radiative corrections. To 
highlight this feature we study the region 7 = i + w for small v in Fig. [31 We 
recognize the familiar feature for n = where the double maxima are replaced 
by a single one a.t v = when the matching to coUinear evolution is performed, 
this removes unphysical oscillations in the gluon Green's function [TU]. We also 
see that all the n ^ Q contributions have a single maximum at that point and 
the resummation only shifts their asymptotic intercepts by a very small amount. 

We would like to point out that the only conformal spin with positive asymp- 
totic intercept is n = 0. This implies that the dependence on the azimuthal 
angles in the Green's function decreases with energy. This property remains 
when the value of the coupling increases as can be seen in Fig. |4l 

After having introduced the resummed kernel used in our calculations we 
are now ready to compare with the experimental results at the Tevatron and 
make predictions for other colliders. This is done in Section [D 



4 Phenomenology 

When dealing with Mueller-Navelet jets we are interested in hadron-hadron col- 
lisions where two jets are tagged in the very forward and very backward regions 
with similar semihard transverse momenta, , such that s ^qcd- 
For large rapidity differences between these two jets logarithms of the form 
(ofg log (s/p^)) should be resummed using the BFKL equation. Mueller and 
Navelet proposed this process in Ref. [3j as ideal to apply the BFKL formalism 
and predicted a power-like rise for the cross section. However, to realize this 
growth as a manifestation of multi-Regge kinematics is very difficult since it is 
drastically damped by the behavior of the parton distribution functions (PDFs) 
for X — > 1. A possible way out is to fix the PDFs and to vary the center-of-mass 
energy of the hadron collider itself, and thereby vary the rapidity difference, Y, 
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Figure 2: Eigenvalues of the scale invariant sector of the BFKL kernel as a 
function of 7 for different values of the conformal spin. 
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between the two tagged jets. BFKL predicts a behavior of the cross sections 
of the form a ^ exp (a — 1)Y /^/Y with a being the intercept. The D0 collab- 
oration analyzed data taken at the Tevatron pp-collider from two periods of 
measurement at different energies yfs = 630 and 1800 GeV. From these they 
extracted an intercept of 1.65 ± .07 [14]. This rise is even faster than that pre- 
dicted in the LO BFKL calculation which for the kinematics relevant in the 
D0 experiment yields an approximated value of 1.45. It has been argued [15] 
that the exact experimental and theoretical definitions of the cross sections dis- 
agreed making the interpretation of the results cumbersome, and the fact that 
the experimental determination of the intercept is based on just two data points 
leaves room for other possible explanations. 

In this work we focus on a more exclusive observable, namely, the azimuthal 
angle decorrelation between the jets. We would like to recall that the Mueller- 
Navelet jets lie at the interface of coUinear factorization and BFKL dynamics. 
The partons emitted from the hadrons carry large longitudinal momentum frac- 
tions and, after scattering off each other, they produce the two tagged jets. 
Because of the large transverse momentum of these jets, the partons are hard 
and obey coUinear factorization. In particular, their scale dependence is gov- 
erned by the DGLAP evolution equations. Between the jets, on the other hand, 
we require a large rapidity difference, a configuration largely dominated by 
multi-Regge kinematics. Therefore, the hadronic cross section factorizes into 
two conventional coUinear parton distributions convoluted with the partonic 
cross section, described within the BFKL context. With respect to the partonic 
cross section, the incoming partons, consequently, are considered to be on-shell 
and coUinear to the incident hadrons. 

For the angular correlation theoretical predictions from LO BFKL were first 
obtained in Refs. [THIII^, improvements due to the running of the coupling and 
proper treatment of the kinematics have been implemented in Refs. USJ '19j. A 
first step towards an analytic NLO description has been made in Refs. |2ll20] on 
which our work builds up. Ten years ago the D0 collaboration at the Tevatron 
measured the azimuthal decorrelation between Mueller-Navelet jets [21]. At 
that time only the LO BFKL equation was available and predictions based on 
it failed to describe the data since it estimates too much decorrelation. Mean- 
while, an exact fixed NLO (af) Monte Carlo calculation using the program 
JETRAD f22l underestimated the decorrelation. In contrast, the Monte Carlo 
program HERWIG [23] was in perfect agreement with the data. 

The main target of our present work is try to improve the prediction for 
this observable using the BFKL resummation introducing NLO effects in the 
kernel. We now show the effect of these new terms for the different observables 
related to the azimuthal angle dependence. As mentioned before, the convenient 
choice of the rapidity variable Y = lns/piP2 turns the convolution with the 
effective parton distributions into a simple global factor which cancels whenever 
we study ratios of cross sections or coefficients C„. For such observables, the 
hadronic level calculation does therefore not differ from the partonic one in 
this approximation. We start by showing in Fig. [5] the Tevatron data for the 
average of the azimuthal angle between the two tagged Mueller-Navelet jets, 



12 



(coS(/)) — Ci/Co and (cos20) = C2/C0, and compare them with our resummed 
prediction developed in the previous section using Eq. p5)) . which evaluates the 
angular mean values in terms of the coefficients C„ in Eq. ((S]). For comparison 
we also show the LO and standard NLO BFKL results without any further 
resummation of higher order terms. As a general trend a decrease of the amount 
of correlation as Y gets larger is obtained, and it can be seen that the NLO 
corrections to the BFKL kernel change the LO results significantly. For the 
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Figure 5: (cos0) = Ci/Cq and (cos 20) — C2/C0 at a pp collider with y/s = 
1.8 TeV for BFKL at LO (sohd) and NLO (dashed). The results from the 
resummation presented in the text are shown as well (dash-dotted). Tevatron 
data points are taken from Ref. 24J. 

particular cuts at the Tevatron, where the transverse momentum for one jet is 
20 GeV and for the other 50 GeV, it turns out that the NLO calculation in the 
MS-scheme provides the best fit to the data. However, this prediction is very 
instable under a change of rcnormalization scheme and we cannot trust it. A 
first hint of this point is that if wc change from MS to GB scheme we notice 
that the NLO result varies more than the LO one. Meanwhile, the resummed 
prediction does not change. This can be clearly seen in Fig. [6] where we have 
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calculated (cost/)) in both renormalization schemes. 

It is important to indicate that the convergence of our observables is poor 
whenever the coefhcient associated to zero conformal spin, Co, is involved. If we 
eliminate this coefficient by calculating the ratios defined in Eq. then the 
only dependence is on the higher n's and the predictions are very stable under 
the introduction of higher order corrections. This is illustrated in Fig. [7] where 
we can observe that the predictions at LO, NLO and with a resummed kernel 
for = ^ are very similar. 

<cos0> Ci 

We have also studied the full angular dependence by investigating the differ- 
ential angular distribution as given in Eq. P71) . The D0 collaboration published 
their measurement of this normalized angular distribution for different rapid- 
ity differences in Ref. [21] . In Fig. [8] we compare this measurement with the 
predictions obtained in our approach using a LO, NLO, and resummed BFKL 
kernel. This comparison is very useful to further justify the need of a collinear 
resummation to all orders. The NLO result here presented is again in the 
MS-scheme, when we switch to the GB-scheme the plot completely changes 
becoming even negative as we approach 4> i'"'- This is not the case in the 
collinearly improved calculation. We can also see that the fit to the data in 
the resummed case is much better than at LO and we have checked that the 
analysis of x^/n.d.f. for the resummed kernel improves for larger rapidities. 
Although for the low rapidities measured at the Tevatron our calculation is not 
close to the data, the fact that the shape of the distribution is the correct one 
is very reassuring. It would be very interesting to have measurements of this 
observable at the future LHC at CERN where a much larger center ~of-mass 
energy is accessible to investigate if a BFKL-based analysis fits the data better 
for larger rapidity differences. This available rapidity range is restricted rather 
by the geometry of the detector than by the energy of the colliding particles 
and by placing calorimeters far enough in the forward and backward regions it 
would be possible to reach about 10 — 12 units of Y difference which would be 
very useful to gauge the importance of multi-Regge-kinematics. 

If large values of Y were accessible in the data it would be very interesting to 
propose other observables where BFKL effects should be visible. As an example, 
with ^12 units of rapidity between the most forward and most backward jets 
wc could tag another jet in the central region of the detector with the condition 
that the three jets had similar transverse momenta of ~ 10 GeV. The two regions 
at Y ~ 6 from the central jet would be enough for BFKL evolution. Studies of 
the growth with energy of this configuration together with double differential 
cross sections in the relative azimuthal angles between the three jets would help 
disentangle the underlying BFKL dynamics. 

We present numerical estimates for the ratios Cm/C„ for a broader range of 
rapidity as predictions for the LHC in Fig. [9l For the differential cross section 
we also provide results at large Y in Fig. [101 Our calculation is not exact and 
we partially estimated the uncertainty associated to the running of the coupling 
and to the sector of the NLO Mueller-Navelet jet vertex originating from the 
splitting functions, which can be easily read off from Refs. [5]. It turns out 
that the effect on the overall normalization can be large, as it has been shown 
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Figure 6: The same plots as in Fig. [5] for (cos (j)) comparing the MS renormal- 
ization scheme (sohd) with the GB scheme (dashed). The plots correspond to 
LO (top), NLO (middle) and coUinearly resummed (bottom) BFKL kernels. 
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1 2 3 4 5 6 

Fignre 7: ^j^^ = §7 as obtained from Fig. E] with LO (solid), NLO (dashed) 
and coUinearly resummed (dash-dotted) BFKL kernels. 

in Ref. [25] , but the influence in the ratios that we consider here is only of the 
order of a few percent. We also estimate the uncertainty in our choice of rapidity 
variable Y by varying our choice of Regge scale sq = piP2 by a factor of 2. We 
also varied the renormalization scale /i by the same factor, and represented 
the uncertainty stemming from these two sources by gray bands in Fig. [9] and 
Fig.Hni 

5 Conclusions 

We have presented a detailed analytic study of the effects of NLO corrections 
to the BFKL kernel in the description of azimuthal angle decorrelations for 
Mueller-Navelet jets in hadron colliders. From the technical point of view we 
found that the intercepts for conformal spins other than zero have good conver- 
gence properties and are not largely modified when a coUincar resummation is 
performed. The zero conformal spin component does need of this resummation 
to get stable results. It is possible to reduce the uncertainties present in our 
analytic study by using Monte Carlo techniques [26, and work is in progress in 
this direction p7] . 

From the phenomenological side we have performed comparisons to the only 
available data related to this observable, which were extracted at the Tevatron 
many years ago. As the rapidity differences between the two tagged jets in the 
D0 experiment were not very large it is natural to think that a fix order calcula- 
tion would do a better job in fitting the data than a BFKL resummation. Indeed 
we find that our results improve with respect to the LO BFKL predictions but 
still show too much azimuthal angle decorrelation. We would like to encourage 
the experimental study of this observable at the LHC at CERN with rapid- 
ity differences quite larger than the presently available in the literature. This 
would be very useful to investigate the importance of BFKL effects in multijet 
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Figure 8: \a.&pp collider at ^/s~\.d> TcV using a LO (stars), NLO (squares) 

and resummed (triangles) BFKL kernel. Plots are shown for Y = 1 (top), Y = 
3 (middle) and Y = 5 (bottom) . Tevatron data points are taken from Ref. [21] . 
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Figure 9: Different ratios of the coefficients C„ obtained using a coUinearly 
resummed BFKL kernel. The gray band reflects the uncertainty in sq and in 
the renormalization scale /x. 
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production in a hadronic environment. If the available rapidity range is large 
enough there is a good opportunity to study other jet topologies which could be 
dominated by multi-Regge kinematics. For these studies the NLO production 
vertex and techniques developed in Ref. [71 [SH] will be very useful. 
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